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A thorough mapping between the hydrodynamics of a two-dimensional Bose-Einstein condensate 
and the nonrelativistic classical electrodynamics of a charged material medium is proposed. This is 
shown to provide a very useful frame to discuss several features of vortex dynamics. Two important 
local conservation theorems of energy and angular momentum are derived and further applications 
are summarized. 



I. INTRODUCTION 

Popov pj was the first to explain how the hydrodynamics of vortices in two-dimensional (2D) homogeneous super- 
fluids can be mapped onto (24-l)D relativistic electrodynamics, with vortices playing the role of charged particles and 
phonons the role of photons. The dynamics of superfluid films was later investigated by exploiting the analogy with 
the dynamics of a 2D plasma Q. The mapping onto (2+l)D electrodynamics was utilized in Ref. 3] to explore the 
theory of dynamical vortices in superfluid fllms, deriving a frequency dependent vortex mass. Such a mapping was 
also employed to study the phonon radiation by a vortex, which exercises circular motion under the influence of an 
external force in a 2D homogeneous superfluid Q- More recently, Fedichev et al. [sj, using the analogy to electron- 
positron pair creation in quantum electrodynamics, investigated the process of vortex-antivortex pair creation in a 
supersonically expanding and contracting quasi-2D Bose-Einstein condensate (BEC) at zero temperature. However, 
in spite of the current interest and widespread importance of the phenomenon of Bose-Einstein condensation of dilute 
gases in traps, the analogy of 2D BECs with electrodynamics remains almost unexplored. 

We propose in this paper a thorough mapping of the dynamics of a 2D BEC onto the electrodynamics of a 2D 
material medium. After reviewing the main features of a 2D BEC in Sec. [TTl we present such a mapping in Sec lIIII 
Next, this is applied in Sees. IIVI and IVl to deriving local conservation theorems of energy and angular momentum, 
respectively. Finally, further applications are summarized in Sec. IVII 



II. TWO-DIMENSIONAL BEC 

Our starting point is the effective 2D Gross-Pitaevskii equation Q 

z;.|ci>(r,0= (^-^ + yext(r)+5|$(r,0p)$(r,t), (2.1) 

where 14xt(r) denotes the trapping potential and g corresponds to the effective 2D coupling constant between the 
atoms of mass m. The complex order parameter $ may be written in terms of a modulus and a phase, as follows: 



$(r,i) = v/'^(r,t)exp[i5(r,t)], (2.2) 
where n(r, t) denotes the particle density and the gradient of the phase S yields the velocity field, 

v(r,t) = -V5(r,i). (2.3) 

m 
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Such a field is irrotational except at points Vj, where the phase presents a singularity corresponding to a quantized 
vortex, i.e. 

V X V = ^ (5(r - Tj) z, (2.4) 

j 

Kj being the circulation of the velocity field along a path encircling a single vortex located at . 

The Gross-Pitaevskii equation ()2.ip is equivalent to the following two coupled equations for the density and the 
velocity field 0: 

S = -V-(vn) (2.5) 



Assuming that the repulsive interaction among atoms is strong enough, the density profiles are smooth and one can 
safely neglect the kinetic-pressure term, proportional to It? , in the last equation, which then takes the form 



dt \ m m 2 



(2.7) 



This result corresponds to the equation of potential flow for a fluid whose pressure and density are related by the 
equation of state 

p = gn^/2. (2.8) 



The approximation ()2.7p . however, overlooks the vortex dynamics represented by the time derivative of eq. ()2.4p . 
which may be explicitly taken into account by adding a Magnus force term Q: 



^ ' 3 

The equation of continuity (j2.5p and the momentum equation (j2.9p represent a system of hydrodynamic-type equations, 
which after linearization, yields sound waves propagating at the local sound velocity 




(2.10) 



Having explicitly introduced the vortex coordinates Tj into the condensate equations of motion, one must consider 
additional equations given by vanishing Magnus forces acting upon each vortex: 

n'(r,)mK,zx [r,-v'(r,)] =0, (2.11) 

where the prime on particle density and velocity field denotes that the self-contribution of the corresponding vortex 
must be ignored. Note that the above equation simply states that the vortex core will move with the background 
superfluid velocity. However, various effects causing possible departures from this behavior have been pointed out 



III. MAPPING ONTO ELECTRODYNAMICS 



The above two-dimensional dynamics can be mapped onto the electrodynamics of a 2D material medium, with 
electromagnetic waves representing sound waves and macroscopic (free) point charges representing vortices. We shall 
utilize in the following the Heaviside-Lorentz system of units 12]. To begin with, we assume a transverse magnetic 
induction proportional to the condensate density 



(3.1) 
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and an electric field, whose expression arises from representing vortices as point charges of negligible mass, leading to 
a vanishing Lorentz force 



E+-xB = 0, (3.2) 

c 



where c — \J g nraax/fn denotes the maximum speed of sound in the condensate. The above expression, rewritten as 
E/B = w/c, embodies an important analogy to quantum electrodynamics. In fact, the Schwinger vacuum breakdown 
[13| is a phenomenon occuring whenever the electric field exceeds the magnetic field, which gives rise to electron- 
positron pair creation. On the other hand, it is well known that a supersonic superflow becomes unstable with respect 
to the spontaneous creation of vortex-antivortex pairs 5] . 

We shall restrict ourselves to slow motion superflows v/c <Si 1, which ensures a mapping onto nonrelativistic classical 
electrodynamics. Then, from eqs. p.ip and p. 21) we have 

E^^nzxv. (3.3) 

c 

Equation (13.11) implies the absence of free magnetic poles, V • B = 0, while the continuity equation (12. 5p becomes 
equivalent to Faraday's law: 

VxE.-if^ (.4) 

The electric displacement field may be obtained by assuming a permittivity given by. 



71,1 



e 



which, according to (13. 6p and (13. 8p . yields 



(3.5) 



c; n 
and therefore 

D = e E = rimax z x v. (3.6) 
Thus, the curl equation (|2.4[) yields Coulomb's law: 

V-D = p = ^g,<5(r-r,), (3.7) 

i 

where a negative charge corresponds to a positive vorticity and viceversa. 



Qj = - \/m Umax tj- (3.8) 
Particularly, the velocity field of the vortex located at Vj reads 



2 X (3.9) 



"^-2. |r-r,P' (3-10) 



that is, the correct expression for the field of a point charge qj in two dimensions. 
Finally, the momentum equation (12. 9p yields Ampere-Maxwell law, 

VxH-J/c=i^, (3.11) 
c at 

where the current density is given by 

J = ^g,J(r-r,)r„ (3.12) 
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which is in agreement with the charge density of the right-hand side of eq. (|3.7|) . The magnetic field is H = B — M 
with a magnetization given by 

that is, a permanent magnet perturbed by a constitutive second-order correction in the displacement field. 

The above Maxwell equations must be suplemented with vanishing Lorentz forces acting upon each point charge 
(cf. eqs. (|21T|) and dSj])), 



1j 



E'(r,) + ^xB'(r,) 



0, (3.14) 



where the primed fields denote that they are excluding the self-contribution of the charge qj . 

IV. CONSERVATION OF ENERGY 

Now we shall analyze Poynting's theorem 12]. The total rate of doing work by the fields is J • E per unit area. 
This power must be balanced by a corresponding rate of decrease of energy in the electromagnetic field. In order to 
exhibit this conservation law explicitly, we first use the Ampere-Maxwell law p.ll|) to eliminate J, 

aD , , 

J-E = cE-V X H-E- — . 4.1) 

dt ^ ' 

Now, employing the vector identity E-VxH = H- VxE — V-(ExH) and Faraday's law p.4|) . the right side of 
(|4.ip becomes 

- --(B • H - SV2) - cV • (E X H), (4.2) 

where the last line stems from p.lSp and p.6|) . Here we may recognize the differential form of Poynting's theorem 

^-t-V.S = -J.E (4.3) 
at 



with 



« = B . H - SV2 = n ( ^ + + ^ ) (4.4) 



the condensate energy density, and the vector 



S = cExH = — nv, (4.5) 
on 

representing energy fiow, which may be called the BEC Poynting's vector. 
Taking into account (|3.12p . the right side of eq. (14. 3p reads 

- J • E = - ^ <5(r - r,) E'(r,) • r„ (4.6) 

where we have disregarded the self-field contribution to the work. Two important conclusions may be drawn from the 
above expression. First, by comparing eqs. (j3.14p and (j4.6p . we may see that Poynting's theorem for the condensate 
simply reads 

^+V-S-0, (4.7) 
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which means that energy can only change through the flow of Poynting's vector. The second important expression 
arises as follows. From (|4.3p and (|4.6p we have 

du 

— = -q,S{r-r,)-E'{r,), (4.8) 

which integrated over the condensate yields, 

1 dU 

being U the condensate energy. Then, inserting the above expression into eq. (|3.14[) we obtain 

1 dU . 



E'(r,) = -- — , (4.9) 



mKjnQ{rj) dvj 



(4.10) 



where the condensate density in absence of the vortex, n'(rj), was approximated by the density of the ground state 
no(rj). Expression (j4.10p has proven very useful for numerically evaluating the vortex velocity in terms of the 
condensate energy gradient, yieldin g ac curate results even in cases where the vortex velocity differs appreciably from 
the background superfluid velocity . 

V. CONSERVATION OF ANGULAR MOMENTUM 

In this section we shall utilize the electromagnetic analogy to derive a local conservation law of angular momentum 
for a BEC in an axisymmetric trap. First, using the Maxwell equations, we may write the electromagnetic force per 
unit area as (T2j . 

1 Id 
pE+-JxB = E- VD-Dx(VxE)-Bx(VxH)-- — (DxB). (5.1) 

c c at 

Then, taking into account the axial symmetry of the trapping potential, the torque per unit area may be written 



rx(/9E+ijxB) = i\x 



- y 

_i|[rx(DxB)]. (5.2) 
c ot 

On the other hand, from Newton's second law, such a torque corresponds to the rate of change of the density 
of mechanical angular momentum £mcch- Therefore, eq. (|5.2[) may be rewritten as the differential form of the 
conservation law of angular momentum: 

I- {Crncch + Acid) + V • X = 0, (5.3) 

ot 

where the field angular-momentum density is 

'Cficid = r X (D X B)/c = rtr x (mv), (5.4) 
i.e., the BEC angular-momentum density. The flux of angular momentum in (|5.3p is described by the pseudo tensor 

m 

M=Tx r, (5.5) 
where the tensor T has nonvanishing components only in the x-y subspace 

= —mnViVj — 6ij p, (5-6) 
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with the pressure p given by eq. (j2.8|) . Here we may recognize —T as the BEC momentum flux density tensor. 

Next we show that the torque (|5.2p vanishes. To this aim, we may rewrite the force given by the left side of (|5.ip 
using eqs. (|3.7|) and (|3.12|) . which yields 



E'(r,) + ^xB'(r,) 



(5.7) 



where we have excluded self-force contributions. But, according to (|3.14p . the above expression between square 
brackets must vanish, which yields a vanishing torque. Hence, the BEC conservation law of angular momentum 
finally reads 

d 

— (nmr X v) + V • X = 0, (5.8) 
which means that angular momentum can only change through the flow of the tensor A4. 

VI. FURTHER APPLICATIONS 



The following applications of the present formalism are treated in Ref. 



• The velocity field induced by an off-centered vortex, which is analyzed from the viewpoint of the equivalent 
electromagnetic problem. 

• Two intriguing aspects of the vortex dynamics: 

— The recently found 'core effect' on the vortex velocity which is interpreted within the electromagnetic 
picture. 

— The highly controversial question of the vortex inertia, which is discussed from the viewpoint of the elec- 
trodynamic analogy. 
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